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Abstract. Multiple sets of synthetic spectra of OB- 
binary stars are used to test the suitability of disentan- 
gling for deriving accurate spectroscopic orbits. Given a 
set of spectra with broad phase coverage and sufficient to- 
tal integration time (almost independent of the number 
of spectra), it appears that disentangling yields accurate 
and reliable semi-amplitudes of the spectroscopic orbits 
(other parameters being fixed). Methods for estimating 
the uncertainties on the derived semi-amplitudes are in- 
vestigated. 
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1. Introduction 

The disentangling technique was first presented by Simon 
& Sturm (1994, hereafter SS94) and applied to the early Re- 
type binary, V453 Cyg. It has subsequently been used on 
the O-type systems DH Cep (Sturm & Simon 1994) and 
Y Cyg (Simon et al. 1994). The original motivation for 
the method was that it enabled the separation of closely 
blended spectra, which were then suitable for quantita- 
tive analysis. Accurate separation of the spectra requires 
accurate parameters for the spectroscopic orbit. Those pa- 
rameters not fixed by other information can, in principal, 
be determined by finding the values of these free parame- 
ters that minimize the residual of the fit to the observed 
spectra when the separated spectra are recombined. 

In this research note we summarize the results of 
Hynes (1996), in which multiple sets of synthetic spectra 
were used to investigate whether disentangling is a viable 
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method for determining accurate spectroscopic orbits for 
a selection of OB eclipsing binary stars. 

2. The synthetic data 

The synthetic binary star spectra are based on a grid 
of single star spectra kindly provided by Dr C.S. Jeffrey 
(1996). These cover the spectral range 3600-5000 A. Each 
component spectrum was obtained by linear interpolation 
on the single star grid using the best estimates of T e ff 
and logg, followed by convolution with the appropriate 
rotational broadening profile. A Doppler shift was then 
applied to each component based on the adopted spec- 
troscopic orbit and the phase desired and the two shifted 
spectra were added according to luminosity ratio. Pois- 
son noise was applied to the spectra to produce a realistic 
signal-to-noise ratio (S/N) (typically S/N= 100-200). Ten 
data sets (each of 15-20 spectra) were generated for each 
case studied, with phases randomly chosen from a uniform 
distribution excluding eclipses. Parameters of the systems 
discussed here are given in Table |l|. 

For eclipsing systems parameters other than the semi- 
amplitudes of the spectroscopic orbit (Ki and K2) can 
usually be determined from photometric data. We can 
therefore estimate optimum values of K\ and K2 for each 
of the ten synthetic data sets using a robust grid searching 
method. Briefly, the residual, r is calculated for a range 
of K\ and K% (using our own implementation of SS94's 
algorithm), using successively finer grids centered on the 
maximum of the previous grid. The mean difference be- 
tween the ten resulting K\, K% values and the adopted 
values indicates the reliability of the method (i.e. this is 
a test for systematic errors) and the population standard 
deviation of the ten sets measures its accuracy. 

3. Phase distribution 

Existing data sets for binary stars are often concentrated 
around phases near quadrature to minimize problems with 
blending inherent with established techniques for measur- 
ing radial velocities. This is not the best approach to take 
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Table 1. Adopted parameters of the systems studied together 
with mean velocities (and errors on this mean) derived by dis- 
entangling 10 sets each of 20 synthetic spectra. The final col- 
umn shows the population standard deviation of the 10 sets 
and is a measure of the intrinsic scatter of the method. 
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with the disentangling procedure. It is clear from the way 
the equations are set up that the system is best condi- 
tioned if the phases cover the full velocity range. This is 
confirmed by numerical tests (Hynes 1996). Thus it ap- 
pears that, in contrast to other techniques, the disentan- 
gling process will work best when given as broad a phase 
coverage as possible. This will often be to the observer's 
advantage. 

4. Integration time 

We next consider the effects of varying the number of spec- 
tra obtained and the duration of each exposure. Both of 
these can be thought of as varying the total integration 
time on the target. The combinations chosen are listed in 
Table |^, together with the resulting uncertainties in de- 
duced parameters, again, for star C. They are plotted in 
Fig. |l| as a function of total integration time for the pri- 
mary velocity; similar curves are obtained for K2- Two 
power law best fits are plotted; one for variations in num- 
ber of spectra and the other for variations in exposure 
time. 

It is clear from this figure that there is no signifi- 
cant difference between the two curves i.e. it matters little 
whether we change the number of spectra obtained or the 
exposure time of each; it is the total exposure time that 
matters. 



5. Accuracy and reliability of the derived orbits 

In Table [l] we show the adopted K\, K2 values and the 
mean deduced values derived from 10 sets of synthetic 



Table 2. Effect of varying the number of spectra and the in- 
tegration time for Star C. 
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Fig. 1. The effects of varying the exposure time and num- 
ber of spectra on the uncertainty in the deduced primary 
semi-amplitude, K\, Star C. Similar results are obtained for 



spectra. All the systems studied, with the exception of 
star F, show no systematic difference between the adopted 
and deduced K\, K% values i.e. disentangling yields reli- 
able orbital parameters for systems with vsini K. The 
accuracy of disentangling is estimated from the standard 
deviation of K\, K2 between the 10 sets (Table [l], final 
column). Accuracies ^ 1% are found with the exception 
of star F and star D, presumably due to the high rota- 
tional velocity and extreme luminosity ratio, respectively. 
We also looked at the effect of rectification errors ~ 5% 
on the accuracy and reliability of deduced K\ , K2 values 
but found no significant effect. 

6. Estimation of errors 

In practice, given a single real data set, we require some 
method to estimate the accuracy of the K\ , K2 values de- 
rived from disentangling. In this section we test two pro- 
posed methods. 



R.I. Hynes & P.F.L. Maxted: Disentangling as a method of deriving spectroscopic orbits 



3 



6.1. Curvature analysis 

As outlined by SS94 and elaborated by Sturm (1994), we 
can estimate the uncertainties in the derived parameters 
by converting r to x 2 using the mean S /N of the observed 
spectra. The curvature of the x 2 surface (as a function of 
the unknown parameters) then yields the uncertainties in 
the free parameters. The most general method of obtain- 
ing confidence intervals from this surface is to construct 
a contour of constant x 2 ; defined by \ 2 — Xmm + 
and 'project this onto the K\ and axes. For the lcr, 
two parameter case A\ 2 = 2.3. This assumes a normal 
distribution of errors which is a good approximation for 
Poisson noise in the high signal-to-noise limit. 

We determined this contour for star E using 10 sets 
of 20 spectra. The confidence limits on K\ and K2 are 
shown in Table [|. While the deduced uncertainties for K\ 
are satisfactory, those for K2 are somewhat low. This is 
not unique to this data set - it is a common problem. 



6.2. Multiple spectra method 

This study was undertaken in preparation for observations 
using an echelle spectrograph. We therefore decided to in- 
vestigate whether the standard deviation of the results 
obtain with different orders of the cchellogram would pro- 
duce reliable uncertainties. We decided upon a simplistic 
approach, with each order weighted equally. This will de- 
grade the results somewhat, but has the advantage that it 
does not involve estimating the error from a single order 
and is therefore independent of methods using the curva- 
ture matrix. The spectral range of our synthetic spectra 
was equivalent to only 4 echelle orders. For the specific 
case of star E, Table [| shows that the results obtained 
by dividing spectra into four sections, each analyzed in- 
dependently (note that this used a different data run to 
that in Table |^, so the choice of random phases and the 
noise will be different). On average, this is an unbiased 
estimate of uncertainties but the low number of "indepen- 
dent trials" results in poor estimates of the uncertainty in 
an individual case. 



Table 3. la confidence limits from the projection of the 
X 2 = Xmin + 2-3 contour for star E. The population standard 
deviations between the data sets in this table are 0.25 for A"i 
and 0.51 for K%. 
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Table 4. Results of disentangling ten synthetic data sets in 
four separate segments. The population standard deviations 
between the data sets in this table are 0.4 for K\ and 0.7 for 
K2. Those obtained by disentangling the whole spectra were 
0.25 for Ki and 0.48 for K 2 . 
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This is not quite the approach taken by SS94 who 
instead measure the curvature of the surface at the 
minimum, represented by the curvature matrix, T. The 
confidence intervals, 5Ki, are then given by SKi = 
iv/AxVO^" 1 )" (Sturm 1994, Press et al. 1992). We de- 
termine the curvature numerically using the data from the 
final grid for which points are separated by 0.1 kms -1 , 
which is similar to the expected uncertainty i.e. a suitable 
value for determination of the curvature. Provided that 
the surface is smooth, the resulting uncertainties are com- 
parable to those obtained by projection of the confidence 
region. 



7. Discussion 

Our experience with disentangling has led us to be cau- 
tious about uncertainties determined from the curvature 
of x 2 space near its minimum. These rely on several as- 
sumptions which may not be valid in the case of disen- 
tangling. For example, it is not clear that the disentan- 
gling process is strictly equivalent to model-fitting by least 
squares minimization and that we can treat it as a two- 
parameter (i.e. K\ and K2) problem. The disentangling 
procedure simultaneously fits both separated component 
spectra and the radial- velocity amplitudes; in effect each 
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Fig. 2. Two grids resulting from disentangling synthetic spectra of LZ Cen differing only in their random elements. 



pixel of each separated spectrum is a parameter to be 
fitted. In applying curvature analysis as above we are as- 
suming that the separated spectra are simply "nuisance 
parameters" and can be ignored in the error analysis. This 
may seem a reasonable assumption but the fact that cur- 
vature analysis consiste ntly gets the relative errors on K\ 
and K 2 wrong (Section |6j| ) indicates that the role of the 
separated spectra may be more subtle. 

Statistical subtleties aside, a more pragmatic reason 
for caution is shown in Fig. |[ Two grids are shown for 
two synthetic data sets of the same star that differ only 
in their noise and phase distribution (which are similar). 
One grid shows a "ridge" near the minimum, a feature 
that was commonly seen in these grids. The uncertainties 
derived from the curvature of such a surface are clearly 
unreliable. In this case, constructing a surface of constant 
X 2 would be more appropriate. In surfaces derived from 
real spectra, even more complex structure can sometimes 
be seen, for example multiple maxima or a series of parallel 
ridges. Such effects will further complicate the analysis. 

The multiple spectra method appears to give unbiased 
estimates of the uncertainties but requires more indepen- 
dent spectral sections showing useful spectral lines than 
are seen in early- type stars. It might well be applicable to 
late-type stars where even small spectral regions contain 
a large number of features. 

In attempting to analyze the accuracy of the proce- 
dure with synthetic data sets we must consider how well 
we reproduce the characteristics of real data. One issue 
is that the phase sampling achieved is likely to be less 
than ideal. We have attempted to mimic this by choos- 
ing random phases with no extra constraints (for example 
we do not impose any restriction on how close in phase 
two observations can be - real observations spread over 
several nights could give duplicate phasing). The sharp 
ridges commented on above may be one manifestation of 
poor sampling. 



Another way in which real data will differ is that the 
noise distribution is unlikely to be the idealized uncor- 
related Poisson noise assumed in this work, e.g. interpo- 
lating spectra onto a logarithmic grid is known to intro- 
duce a short scale autocorrelation into the noise, which can 
lead to a "rippling" on the residual surfaces. Other effects, 
such as cosmic ray events and wavelength calibration er- 
rors may also lead to noticable distortions of the surfaces 
which must be identified and accounted for. Numerical 
simulations of the type discussed in this paper will be a 
valuable tool in investigating these issues. For real data 
however, we will not have the luxury of obtaining multi- 
ple duplicate data sets and comparing them. The role of 
simulated multiple data sets is then to give us insight into 
the reliability of methods such as curvature analysis. 
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